Method for direct observation of coherent quantum oscillations in a superconducting 

phase qubit 
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Time-domain observations of coherent oscillations between quantum states in mesoscopic super- 
conducting systems were so far restricted to restoring the time-dependent probability distribution 
from the readout statistics. We propose a new method for direct observation of Rabi oscillations in 
a phase qubit. The external source, typically in GHz range, induces transitions between the qubit 
levels. The resulting Rabi oscillations of supercurrent in the qubit loop induce the voltage oscilla- 
tions across the coil of a high quality resonant tank circuit, inductively coupled to the phase qubit. 
It is the presence of these voltage oscillations in the detected signal which reveals the existence of 
Rabi oscillations in the qubit. Detailed calculation for zero and non-zero temperature are made 
for the case of persistent current qubit. According to the estimates for decoherence and relaxation 
times, the effect can be detected using conventional rf circuitry, with Rabi frequency in MHz range. 

PACS numbers: 03.65.Ta, 73.23.Ra 



I. INTRODUCTION 

As is known the persistent current qubit (phase qubit) 
is one of the candidates as-a-key element of a scalable solid 
state quantum processor .Bel The basic dynamic manifes- 
tations of a quantum nature of the qubit are macroscopic 
quantum coherent (MQC) oscillations (Rabi oscillations) 
between its two basis states, which are differed by the 
direction of macroscopic current in the qubit loop. 

Up till now the Rabi oscillations in time.domaini@ or 
as a function of the perturbation poweitm have been 
detected indirectly through the statistics of switching 
events (e.g. escapes into continuum). In either case the 
probabilitiy P(t), or P(E), was to be obtained and ana- 
lyzed to detect the oscillations. 

More attractive in the long run is a direct detection of 
MQC oscillations through a weak continuous measure- 
ment of a classical variable, which would implicitly in- 
corporate the statistics of quantum switching events, not 
destroying-iii. the same time the quantum coherence of 
the qubit Jffl 

In this paper we describe the approach which allows 
a direct measuring of MQC oscillations of macroscopic 
current flowing in a loop of a phase qubit. This qubit 
variety has the advantage of larger tolerance to external 
noise, especially to dangerous random background charge 
fluctuations.Ell To be specific, we will use the example 
of three-junction .small- inductance phase qubit (persis- 
tent current qubittl) where level anticrossing was already 
observed. EJ 

In our method a resonant tank circuit with known in- 
ductance Lt, capacitance Ct and quality factor Qt is 
coupled with a target Josephson circuit through the mu- 
tual inductance M (Fig. |l|). The method was successful^ 
applied to a three-junction qubit in classical regime,t3 
when the hysteretic dependence of ground-state energy 
on the external magnetic flux was reconstructed in ac- 



FIG. 1: Phase qubit coupled to a tank circuit. 



cordance to the predictions of Ref. g. 

The phase qubit is biased by external magnetic flux 
$(t). Assuming small qubit self- inductance, we neglect 
the shielding current. Therefore the flux through the 
qubit loop is 



$(t) = <$> x + * oc (t), 



(1) 



where <fr x is a time independent external flux, $ ac (t) is 
a monochromatic high frequency flux from the external 
source. 

If time-dependent external flux is applied to the qubit, 
the latter will be in a time-dependent superposition of 
states |0) and |1). If the frequency of external flux is in 
resonance with the interlevel spacing of the qubit, the 
average current in the qubit loop will oscillate with the 
frequency f2 which depends on the amplitude of external 
flux and can be made much smaller than the frequency 
of external signal. These smaller oscillations which are 
called Rabi oscillations can be detected with the aid of a 
high quality tank circuit coupled inductively to the qubit 
loop. 

Below we find the expression for the current in a qubit 
loop in presence of electromagnetic resonant excitation 
for zero and non zero temperatures. We show that for 
typical parameters of three-junction phase qubit the cur- 
rent oscillations in a qubit loop give rise to the voltage 
oscillations across the tank with the amplitude which can 
be at the //V level that makes a direct detection of Rabi 
oscillation possible. 
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II. QUANTUM DYNAMICS OF 3JJ PHASE 
QUBIT 

Quantum dynamics of three-junction phase qubit has 
been studied in detail in Ref. ^. The qubit consists of a 
loop with three Josephson junctions. The loop has very 
small inductance, typically in pH range. It insures ef- 
fective decoupling of qubit from external environment. 
Two Josephson junctions have equal critical current Ic 
and capacitance C, while the critical current and capac- 
itance of a third junction is a little bit smaller, a/p, aC 
where 0.5 < a < 1. If the Josephson coupling energy 
Ej = /c'$o/27r, where $o = h/2e is a flux quantum, is 
much more than the Coulomb energy Eq — e 2 /2C, then 
the phase of a Cooper pair wave function is well defined. 
As was shown in Refs. |l]J^ in the vicinity of $ = $0/2 this 
system has two quantum stable states which macroscop- 
ically differ by the direction of the current circulating in 
a qubit loop. In the absence of high frequency excita- 
tion the quantum properties of the qubit are described 
by Hamiltonian (Eq. (12) in Ref. §): 



Ho 



P 2 
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where P v = 

($o/2tt) 2 2C, Me = ($o/2^) 2 2C(l + 2a); 
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U(f,ip,6) = Ej {2 + a -2cos^cos6» 
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(3) 



Here ip = (tpi + tpz)/2, 9 = (tpi - ¥>2)/2), where <p u p 2 
are gauge-invariant phases of two Josephson junctions 
with equal critical currents. In contrast to Ref. |[ in (^) 
we define the flux bias / = $/&o — \ as a small pa- 
rameter measuring the departure from degeneracy. At 
the degeneracy point / = potential energy (||) shows 
two minima with equal energies So at the points tp = 0, 
9 = ±9 C , where cos6> c = l/2a. The tunneling between 
the minima lifts the degeneracy leading to the energy 
levels E± — Eq ± A where A is a tunneling matrix ele- 
ment between two minima. However, at the degeneracy 
point the current in a qubit loop vanishes, so that it is 
necessary to move a little bit away from this point. In 
order to find the levels in the close vicinity of degeneracy 
point we expand potential energy (^|) near its minima 
taking account for linear terms in / and quadratic terms 
in quantum variables </), 6. With the use of the technique 
described in Ref. ^| we find the following expression for 
the energies of two levels: 



E ± =e ±jE 2 f 2 X 2 (a) + A 2 



(4) 



where we take offset at the degeneracy point, i. e., f=0 
corresponds to $x = &o/2; 



£0 — Ej 




1 + y/2a - 1) 



(5) 
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(6) 



Expression (|4|) differs from corresponding equation in 
Ref. 11 by a factor A(a) which explicitly accounts for the 
dependence of the energies E± on a and Ej /Ec- The 
stationary state wave functions are eigenfunctions of 
Hamiltonian Hq : Hq^± = E±^±. They can be written 
as the superpositions of the wave functions in the flux 
basis, ^Li^r where L, R stand for the left, right well, 
respectively: *± = a±9 L + b±^ R . 



a± 



^/(e + -E ± ) 2 + A 2 



;b± = 



e+-E ± 



y/(e+ - E ± ) 2 + A 2 



(7) 

where e+ = (V L \H \V L ); £_ = ($r\Ho\V r ). For sta- 
tionary states the current circulating in a qubit loop can 
be calculated either as the average of a current operator 
I q = Ic sin(<y9 + 9) over stationary wave functions or as a 
derivative of the energy over the external flux: 

h = ^\h\*±) = ^ = ±Icf^^ (8) 



7T flUjQ 



where tku = E + — E_ . 

Suppose we apply to the qubit the excitation on a fre- 
quency close to a gap frequency lo = hu>o f=a 1 GHz. The 
corresponding perturbation term is then added to Hamil- 
tonian Hq: Hi n t = V((j),9)cos(u)t) where V(<f>,0) in the 
vicinity of the left (right) minimum is as follows: 



2a 2 
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(9) 



(10) 



Throughout the paper we assume / > 0, so that the 
bottom of the left well is higher than the bottom of the 
right well. Accordingly, in Eqs. (||) and (|h]) the upper 
sign refers to right minimum, the lower sign refers to left 
minimum. The quantity f ac in Eq. (^) is the amplitude 
of excitation field in flux units: f ac — & a c/®o- The high 
frequency excitation induces the transitions between two 
levels which result in a superposition state for the wave 
function of the system: = C+(t)*+ + C_(i)*_. 

The coefficients C±(t) are obtained from the solution of 
time dependent Schrodinger equation with proper ini- 
tial conditions for C±(t). We assume that before the 
excitation the system was at the lower energy level: 
C_(t = 0) = l;C+[t = 0) = 0. The corresponding solu- 
tion for C±(t) in the rotating wave approximation is as 
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follows 111 

C-(t) : 
C+(t) : 



"V" 



cos S 2r — z — — — sin s it 



2<1 



.O, 



Ll 1 -ijj- ) sill O/ 



(11) 



the amplitude of time-dependent flux which is coupled 
to qubit from high frequency source f ac = 1 x 10~ 4 . We 
set the flux offset from degeneracy point / = 3.5xl0~ 4 , 
so that huiQ = 2V2A, ujo/2it = 0.85 GHz. For these val- 
ues we obtain from Eq. (12)the Rabi frequency f2 r /27r = 
32 MHz. 



where fL 



VLl. At reso- 



nance (u> = u>q) we get O = Q r . Taking ^ l, ^ r as ground 
state oscillatory wave functions in the left, right well, re- 
spectively, we calculate matrix element (\E' + |y|\f r _) and 
obtain for Rabi frequency: 



= ^-/a C |A(a)|-^- 

n nuio 



(12) 



As is seen from Eq. (pjl) by a proper choice of excitation 
power (f ac in our case) the frequency fl r can be made 
much lower the gap frequency A/h. Now we calculate 
the average current in a superposition state: 

I q = (*(t)|£|*(f)) 

= |C + (t)| 2 (* + |f 9 |vI/ + ) + |C_(t)| 2 (vl/_|f 9 |*_) (13) 



III. THE INTERACTION OF THE PHASE 
QUBIT WITH A TANK CIRCUIT 

The problem of a coupling a quantum object to the 
classical one, which is dissipative in its nature, has no 
unique theoretical solution. A rigorous approach is to 
start from exact Hamiltonian which describes the qubit- 
tank circuit system: 



H — Hq + Ht + Hot + Htb + Hb 



(17) 



where Ho is the qubit Hamiltonian given in Eq. (^]), Ht 
is the tank circuit Hamiltonian 



Ht — 



$ 2 





2Ct 2Lt 



(18) 



where we neglect high frequency loq term. Accounting 
for the expression (^) for the current in stationary states 
we obtain at resonance the following expression for the 
average current flowing in a qubit loop: 



In 



-Icf 



A 2 (a) Ej 

7T TvjJq 



cos 2fL< 



(14) 



At finite temperature the calculation of the current 
in a qubit loop in presence of high frequency excita- 
tion is based on the density matrix equation: ihp(t) = 
[(Hq + H int (t)) , p(t)\ with initial conditions at thermal 
equilibrium: p ++ (0) = p+ + 



p |_(0) = 0, where p e . 



(0) = p+-(0) = 
are equilibrium density 



matrix elements 



P-\ 



P- 



z 



k B T -)- e k B T . 



In the rotating wave approximation the diagonal ele- 
ments of density matrix at resonance (u> = loq) are as 
follows: 

11 ^0 ^ 

p— ( < ) = 2 + 2 t 2fc^r C0S ' 

p++(t) =l-p__(t). 

Now we find the average current at resonance for 
nonzero temperature: 



h = 



= -I, 



*+\Iq\* + )P+ 

Ejf\ 2 (a 



hTuoq 



tanh 



Tiujq 
2k B T 



P-(t) 
cos2Vt r t (16) 



In order to estimate the Rabi frequency fl r we take 
the following qubit parameters: Iq — 400 nA, A/h — 
0.3 GHz, a = 0.8, L — 15 pH, Ej/E c = 100. We take 



where Q and $ are a quantum operators of the charge 
at the capacitor and of magnetic flux trapped by the 
inductor of a tank circuit, respectively. The operators 
obey commutator relation Q] = ih. The interaction 
Hamiltonian between the qubit and the tank, Hot is: 



H 



M 



OT 



(19) 



where M is inductive coupling between qubit and the 
tank; Hb is the Hamiltonian of a thermal bath coupled 
to the tank via interaction Htb — ot$>T, where a is the 
coupling constant between the tank and dissipative envi- 
ronment, r is the dynamic variable of thermal bath Hb- 
The equations of motion for tank circuit variables are 
as follows: 



dt 



$ M f 
Lt L 

~~dt 



aT 



Q_ 

Ct 



(20) 



(21) 



From these two equations we get for the voltage operator 
V = Q/Ct across the tank: 



d 2 V 



luIv = -Mut^t- 
T T dt 



dT 
dt 



(22) 



The averaging of this equation over the bath leads to the 
dissipative equation for the average voltage, V, across 
the tank (see, for example, Ref. nj): 



V 



^V + cotV^-Mco^. 
Qt dt 



(23) 
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where Qt = ojtRtCt 3> 1 is the tank quality factor, 

UlT = 1/V LtCt- ri 

In the spirit of selective quantum evolution approach,ta 
we interpret Eq. |23| as follows. Suppose the qubit is in a 
pure state, o|0) + &|1). The tank voltage is measured 
(in quantum-mechanical sense) at certain times tfc — 
At, 2At, .... Each time the qubit state is also measured, 
since there is correspondence between V(tk) and qubit 
current, which thereby takes value either (0|ig(tfe)|0) 
or (1| I q (tk) |1), with appropriate probability. Averaging 
V(tk) over intervals AT 3> At, which are still small com- 
pared to other characteristic times in the system, would 
yield Eq. ||. 

The quantity to be detected is the oscillating voltage 
across the tank which at resonance is V = V a cos 2f2 r t, 
where the voltage amplitude, V a is: 



Va = M Qicf^-^-n r 

7r nu>o 



(24) 



It is the presence of these voltage oscillations in the de- 
tected signal which reveals the existence of Rabi oscilla- 
tions in the qubit. 

The Eq. ( |23| ) is still quantum equation since the av- 
erage voltage V is a quantum operator of the voltage 
across the tank averaged over the bath degrees of free- 
dom. Therefore, our problem is reduced to the problem 
of measuring a weak external force (in our case, Mdl q /dt) 
by a dissipative oscillator. As is known (see, for exam- 
ple, Ref. [l5]) a classical descriptions of such oscillator 
requires that the quantum fluctuations of the detector 
variable (i. e. V in our case) in the measurement band- 
width be smaller than the amplitude induced in the tank 
coil by external signal, MQxdI q /dt. According to the 
fluctuation-dissipation theorem the quantum fluctuations 
of the voltage V are given by the spectral density: 



S v (uj) = 2ReZ{uj)hcjcoth ( — — ) 

\2fcsT / 



(25) 



where Z~ x (<jS) = (1/Rt + l/icoLr + iojCt) is the 
impedance of a tank circuit. 

Below we take the following parameters of the tank 
circuit: Lt = 50 nH, Q = 1000, inductive coupling to 
qubit, fc 2 = M 2 /LLt = 10~ 4 . We assume the tank is 
tuned to the frequency of oscillating current in the qubit, 
i. e., w t /2tt = 2tt r /2Tr = 64 MHz, hence, C T = 124 pF. 
We take the flux amplitude which is coupled to qubit 
from high frequency source f ac = 1 x 10 -4 . We set the 
flux offset from degeneracy point / = 3.5 x 10~ 4 , so that 
hujo = 2V2A, luq/2-k — 0.85 GHz. For these values we ob- 
tain for the voltage amplitude at resonance V a ~ 0.7 /iV. 

From Eq. ( p5| ) we estimate the voltage fluctuations 
across the tank circuit coil at resonance (u> = lot) ■ 



V Z(ujt) 



(26) 



where B = lot/^Qt is the bandwidth of the tank cir- 
cuit. We perform the calculations for T = 10 mK. For 



the voltage fluctuations we get V n « 10 nV. Thus, we 
see that the voltage fluctuations across the tank coil is 
much smaller than the signal amplitude from the qubit. 
Therefore, we may treat the tank circuit as a classical 
object and the voltage V(t) as the classical variable cou- 
pled through Eq. 23 to the qubit degree of freedom, where 
the current I q in Eq. ^3] is calculated as the average of a 
current quantum operator over the qubit statistical op- 
erator. 



IV. THE EFFECTS OF QUBIT RELAXATION 
AND DECOHERENCE 

The estimations we made above are very promising, 
however, our derivation is made under a strong assump- 
tion: we neglected the decoherence due to interaction of 
the qubit with external environment and with a measur- 
ing device. In fact, the possibility of detection of Rabi 
frequency depends crucially on relaxation, r r , and deco- 
herence, Tcp rates, which lead to the decay of Rabi oscil- 
lations. The detection is in principle possible if period of 
oscillations is small compared to min[l/T v ,l/T r ]. The 
decoherence is caused primarily by coupling of a solid 
state based phase qubit to microscopic degrees of free- 
dom in the solid. Fortunately this intrinsic decoherence 
has been found to be quite weakEj the intrinsic decoher- 
ence times appeared to be on the order of 1 ms which 
is several orders of magnitude more than period of cur- 
rent oscillations we estimated before: 7r/f2 r = 15.6 ns. 
However, the external sources of decoherence are more 
serious. In our method these are the microwave source 
which induces the Rabi oscillations and the tank circuit 
which has to detect them. From this point two structures 
of microwave source has been analyzed: coaxial line that 
is inductively coupled to the qubitO and., on-chip oscil- 
lator based on overdamped DC SQUIDJ13 The analysis 
has shown the relaxation and decoherence rates were on 
the order of 100 fis at 30 mK for coaxial line, and 150 fis 
and 300 /xs for relaxation and decoherence, respectively, 
for on-chip oscillator at 1 GHz. Here we estimate deco- 
herence times that are due to -a. tank circuit using the 
expressions for T r and T v frorrO 



1 



1 



r '- = ^ = o U— J ( UJ o) coth 



2 V TujJi 




2 



2k B T) 



2nr\ 



EjfXY k B T 

hujQ 



(27) 



(28) 



where dimensionless parameter 77 reflects the Ohmic dis- 
sipation. It depends on the strength of noise coupling 
to the qubit. In our subsequent estimations we take 
77 = 5 x 10~ 3 which is relevant for weak damping limit. 
Below we use thfi|-apBroach described in Ref. [2(J in its 
simplified form.EZHiSE 1 ] While it is not mathematically 
rigorous, nevertheless it gives a correct order of magni- 
tude of relaxation times. The quantity J(u>o) is the zero 
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temperature spectral density of the fluctuations of the 
gap energy of the qubit: J(uj) =< Se(u>)5e(uj) > /h 2 . 
The fluctuations Se are due to the flux noise 5f, which 
is supplied to the qubit by a tank circuit. From Eq.(5) 
we get Se = (4E}\ 2 f /fkJ ) Sf, where Sf = M6I/$ . 
The current noise SI in a tank circuit inductance comes 
from two independent parts: Johnson-Nyquist voltage 
noise in a tank circuit resistance with a spectral den- 
sity at T = 0: SVH = 2hujReZ(uS), where Z~ x (u) = 
(1/Rt + 1/iloLt + iluCt) is the impedance of a tank 
circuit, and from a current noise of preamplifier with 
a spectral density Sa- Therefore, we get for J(u)q): 

J(w ) = (4Ej\ 2 fM/h 2 u) <$> ) 2 Si(u ), where Si(lu) = 

/w 2 L 2 is the spectral den- 



S v {w)ReZ{uj) + S a \Z(uj)Y 
sity of a current noise in the tank circuit inductance. 
Since lu » lu t ,Si(u> ) = {2Tilu /Rt + Sa) (^t/wo) 4 - 
For the estimation we take T — 10 mK, Rt = 20 kf2, 
S A = 10~ 26 A 2 /Hz, w t /2tt = 64 MHz, other parame- 
ters being the same as before. We find that the contri- 
bution to the relaxation of a tank circuit noise and of 
preamplifier noise is approximately 0.1 s" 1 and 22 s _1 , 
respectively. Therefore total relaxation is determined by 
preamplifier noise giving relaxation time T r = 45 ms. 
The decoherence rate is dominated by a second term in 
(f28|), which is equal 5.4 x 10 6 s -1 giving decoherence time 
T v = 185 ns, which is approximately ten oscillation pe- 
riods of circulating current. These estimations clearly cussions. 



show the possibility of detection of Rabi frequency in 
MHz range with the aid of conventional rf circuitry. 

Here we did not consider the effect of a tank circuit 
back action on the qubit which can give additional con- 
tribution to the relaxation and decoherence rates. In 
order to reduce the back action effect of a tank circuit 
and to further enhance the MQC signal it may be ad- 
vantageous for detection of Rabi oscillations to use a two 
dimensional array of identical phase qubits coupled to 
a tank circuit. Modern technology allows several thou- 
sands aL weakly coupled phase qubits to be obtained on 
a chip.EJ Preliminary experiments showed that the be- 
havior of macroscopic current in such a system is com- 
pletely analogous to that of longitudinal magnetization 
in NMR: the collective reversal of the persistent currents 
in the qubit loops when sweeping thc^flux bias within a 
degeneracy point has been observedO 

In conclusion, we have shown that Rabi oscillations 
(in MHz range) of the current circulating in a qubit loop 
which are induced by high frequency external source (in 
GHz range) can be detected in MHz range as the voltage 
oscillations in the high quality tuned tank circuit induc- 
tively coupled to the qubit. 
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